Intrinsic gain modulation and adaptive neural coding

In many cases, the computation of a neural system can be reduced to a receptive field, or a set of linear filters, and a thresholding function, or gain curve, which determines the firing probability; this is known as a linear/nonlinear model. In some forms of sensory adaptation, these linear filters and gain curve adjust very rapidly to changes in the variance of a randomly varying driving input. An apparently similar but previously unrelated issue is the observation of gain control by background noise in cortical neurons: the slope of the firing rate vs current (f-I) curve changes with the variance of background random input. Here, we show a direct correspondence between these two observations by relating variance-dependent changes in the gain of f-I curves to characteristics of the changing empirical linear/nonlinear model obtained by sampling. In the case that the underlying system is fixed, we derive relationships relating the change of the gain with respect to both mean and variance with the receptive fields derived from reverse correlation on a white noise stimulus. Using two conductance-based model neurons that display distinct gain modulation properties through a simple change in parameters, we show that coding properties of both these models quantitatively satisfy the predicted relationships. Our results describe how both variance-dependent gain modulation and adaptive neural computation result from intrinsic nonlinearity.Comment: 24 pages, 4 figures, 1 supporting informatio

The what and where of adding channel noise to the Hodgkin-Huxley equations

One of the most celebrated successes in computational biology is the Hodgkin-Huxley framework for modeling electrically active cells. This framework, expressed through a set of differential equations, synthesizes the impact of ionic currents on a cell's voltage -- and the highly nonlinear impact of that voltage back on the currents themselves -- into the rapid push and pull of the action potential. Latter studies confirmed that these cellular dynamics are orchestrated by individual ion channels, whose conformational changes regulate the conductance of each ionic current. Thus, kinetic equations familiar from physical chemistry are the natural setting for describing conductances; for small-to-moderate numbers of channels, these will predict fluctuations in conductances and stochasticity in the resulting action potentials. At first glance, the kinetic equations provide a far more complex (and higher-dimensional) description than the original Hodgkin-Huxley equations. This has prompted more than a decade of efforts to capture channel fluctuations with noise terms added to the Hodgkin-Huxley equations. Many of these approaches, while intuitively appealing, produce quantitative errors when compared to kinetic equations; others, as only very recently demonstrated, are both accurate and relatively simple. We review what works, what doesn't, and why, seeking to build a bridge to well-established results for the deterministic Hodgkin-Huxley equations. As such, we hope that this review will speed emerging studies of how channel noise modulates electrophysiological dynamics and function. We supply user-friendly Matlab simulation code of these stochastic versions of the Hodgkin-Huxley equations on the ModelDB website (accession number 138950) and http://www.amath.washington.edu/~etsb/tutorials.html.Comment: 14 pages, 3 figures, review articl

Complementarity of Spike- and Rate-Based Dynamics of Neural Systems

Relationships between spiking-neuron and rate-based approaches to the dynamics of neural assemblies are explored by analyzing a model system that can be treated by both methods, with the rate-based method further averaged over multiple neurons to give a neural-field approach. The system consists of a chain of neurons, each with simple spiking dynamics that has a known rate-based equivalent. The neurons are linked by propagating activity that is described in terms of a spatial interaction strength with temporal delays that reflect distances between neurons; feedback via a separate delay loop is also included because such loops also exist in real brains. These interactions are described using a spatiotemporal coupling function that can carry either spikes or rates to provide coupling between neurons. Numerical simulation of corresponding spike- and rate-based methods with these compatible couplings then allows direct comparison between the dynamics arising from these approaches. The rate-based dynamics can reproduce two different forms of oscillation that are present in the spike-based model: spiking rates of individual neurons and network-induced modulations of spiking rate that occur if network interactions are sufficiently strong. Depending on conditions either mode of oscillation can dominate the spike-based dynamics and in some situations, particularly when the ratio of the frequencies of these two modes is integer or half-integer, the two can both be present and interact with each other

Why Are Computational Neuroscience and Systems Biology So Separate?

Despite similar computational approaches, there is surprisingly little interaction between the computational neuroscience and the systems biology research communities. In this review I reconstruct the history of the two disciplines and show that this may explain why they grew up apart. The separation is a pity, as both fields can learn quite a bit from each other. Several examples are given, covering sociological, software technical, and methodological aspects. Systems biology is a better organized community which is very effective at sharing resources, while computational neuroscience has more experience in multiscale modeling and the analysis of information processing by biological systems. Finally, I speculate about how the relationship between the two fields may evolve in the near future

25th Annual Computational Neuroscience Meeting: CNS-2016

Abstracts of the 25th Annual Computational Neuroscience Meeting: CNS-2016 Seogwipo City, Jeju-do, South Korea. 2–7 July 201

Implications of cellular models of dopamine neurons for disease

This review addresses the present state of single-cell models of the firing pattern of midbrain dopamine neurons and the insights that can be gained from these models into the underlying mechanisms for diseases such as Parkinson’s, addiction, and schizophrenia. We will explain the analytical technique of separation of time scales and show how it can produce insights into mechanisms using simplified single-compartment models. We also use morphologically realistic multicompartmental models to address spatially heterogeneous aspects of neural signaling and neural metabolism. Separation of time scale analyses are applied to pacemaking, bursting, and depolarization block in dopamine neurons. Differences in subpopulations with respect to metabolic load are addressed using multicom- partmental models

Optogenetic Stimulation Shifts the Excitability of Cerebral Cortex from Type I to Type II: Oscillation Onset and Wave Propagation.

Constant optogenetic stimulation targeting both pyramidal cells and inhibitory interneurons has recently been shown to elicit propagating waves of gamma-band (40-80 Hz) oscillations in the local field potential of non-human primate motor cortex. The oscillations emerge with non-zero frequency and small amplitude-the hallmark of a type II excitable medium-yet they also propagate far beyond the stimulation site in the manner of a type I excitable medium. How can neural tissue exhibit both type I and type II excitability? We investigated the apparent contradiction by modeling the cortex as a Wilson-Cowan neural field in which optogenetic stimulation was represented by an external current source. In the absence of any external current, the model operated as a type I excitable medium that supported propagating waves of gamma oscillations similar to those observed in vivo. Applying an external current to the population of inhibitory neurons transformed the model into a type II excitable medium. The findings suggest that cortical tissue normally operates as a type I excitable medium but it is locally transformed into a type II medium by optogenetic stimulation which predominantly targets inhibitory neurons. The proposed mechanism accounts for the graded emergence of gamma oscillations at the stimulation site while retaining propagating waves of gamma oscillations in the non-stimulated tissue. It also predicts that gamma waves can be emitted on every second cycle of a 100 Hz oscillation. That prediction was subsequently confirmed by re-analysis of the neurophysiological data. The model thus offers a theoretical account of how optogenetic stimulation alters the excitability of cortical neural fields